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Abstract 

Some modern models of neutron star evolution predict that initially large 
magnetic fields rapidly decay down to some saturation value ~ f ew x 10 13 G 



and weaker magnetic fields do not decay significantly (IPons et all 120091 ) . It 
is difficult to check the predictions of this model for initially highly magne- 
tized objects on the time scale of a few million years. We propose to use 
Be/X-ray binaries for this purpose. We apply several methods to estimate 
magnetic fields of neutron stars in these accreting s ystems using the data 



obtained by the RXTE satellite (IGalache et all |2008[ ). Only using the most 



modern approach fo r estimating th e mag netic field strengths of long period 
NSs as proposed by IShakura et all (120111 ) we are able to obtain a field dis- 
tri bution compatible with predictions of the theoretical model of field decay 
of IPons et all feOPOh . 
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1. Introduction 



Neutron stars (NSs) are the final stage of the evolution of massive stars. 
There are different types of y oung isolated NSs: radio pulsars (see a re- 



view m 



Manchester et al.ll2005l). compa c t cent ral X-ray sources in supernova 
remnants ( IPavlov et all 12004] ; Ide Lucal . 120081 ). magnet ars : anomalou s x-ray 



pulsars (AXPs) and so ft gamma-ray repeaters (SGRs) (jMereghettil . 12008 



Rea and Espositol 2011), clo s e-by cooling radioqu iet isolated NSs - the Mag- 



nificent seven (M7) ( IHaberll . 120071 : iKaplanl . 120081 ). rotating radio transients 
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(RRATs) (jKeand . |2010| ). Another large class of observed systems with NSs is 
formed by X-ray binaries. They are subdivided into high mass X-ray binaries 
(HMXB) and low mass X-ra y binaries (LMXB) (see a rey ie w on X-ray bi- 
naries and their evolution in iPostnov and Yungelsonl ( 120061 ); iBhattacharyya 
(|2009h ). 

In order to deeper comprehend the evolution and to follow links be- 
tween different types of NSs we need a better understanding of the basic 
processes related to NS thermal and magneto-rotational evolution. For the 
initially highly magnetized NSs (B > few x 10 13 G), modern models pre- 
dict that the magnetic field decays significantly o n the time scale of sev- 



eral hundred thousand years ( Aguilera et al. . 20081 ) . For magnetars there is 
a strong observational ey idence for magnetic field decay (see, for instance, 
Pons and Geppertl ( 120071 ) and references therein). These sources are very 



young NSs with ages of the order of tens of thousand years or less and with 
magnetic fields 10 14 — 10 15 G. Unfortunately, it is difficult to use magnetars 
for testing models of magnetic field decay on a time scale > 10 6 yrs because 
they become extremely elusive and avoid detection at ages more than about 
a hundred thousand years. 

For "standard" values of the initial magnetic fields (10 11 - 10 13 G) typical 
for radio pulsars the field decay in these models is insignificant: the magnetic 
field diminishes just by a factor of 2. Therefore, the population of radio 
pul sars is ill-suited to test these model predictions. 



Popov et al.l ( 120101 ) carried out a comprehensive population synthesis of 



young NSs (radio pulsars, magnetars, near-by cooling NSs) and showed that 
a model with magnetic field decay can explain the basic properties of all the 
three populations considered using the same initial magnetic field distribu- 
tion. All types of these objects have ages less than ~ million years. It is 
important to find a way to confront model predictions and observations on a 
longer time scale. NSs in some HMXBs can be suitable for testing the evolu- 
tion of magnetic fields on the time scales of several million years. Compact 
objects in these systems have average ages ~ a few million years (their mas- 
sive companions usually have a tot al lifetime less tha n 10-30 million years). 
In the case of wind-fed accretion ( jNegueruelal . 120101 ). the accretion rate is 
not large enough to affect significantly the evolution of the field (as opposed 
to systems with intense disc accretion which produce millisecond pulsars). 
Thus, by estimating magnetic fields of NSs in HMXBs with accretion from 
a wind, it is possible to study the magnetic field evolution on the time scale 
of a few million years. 
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(Liu et al.. 


2006; 


Reig, 


2011) 



known in the Galaxy, 78 in t he Small Magellanic Cloud (SMC), and 16 in 
the Large Magellanic Cloud ( IRaguzova and PopovL l20Q5h . We consider 40 



Be/X-ray binaries in the SMC to estimate magnetic fields of NSs by means 
of different methods. 

The article is structured as follows. In Sec. 2 we describe the model of 
the magnetic field decay that we compare our estimates with. In Sec. 3 we 
present different methods to obtain estimates of magnetic fields of accreting 
NSs. In Sec. 4 we present the observational data used in the paper. Sec. 5 
contains our results, and in Sec. 6 we discuss them. 

2. Magnetic field decay model 

In this section we describe an analytical approximation for t he model of 



field d ecay used in this paper, following the approach suggested by lAguilera et al. 



( 120081 ) and discuss types of objects that can be used for testing this model. 

There are several mechanisms for the magnetic field decay, but the rate of 
field decay is mainly determined by the Ohmic decay and the Hall drift. The 
Ohmic decay involves diffusion of the magnetic field line s with respect to the 



charg ed particles. The characteristic Ohmic time scale is (IGoldreich and Reisenegger 



19921 ): 

T~Ohm — n • UJ 

cr 

Here A is the typical magnetic field length scale and a is the conductivity of 
NS matter, c is the speed of light. 
The Hall drift time scale is: 

_ 4vren e A 2 

THcdl — „ • \ l ) 

cB 

Here n e is electron density and B is magnetic field strength. The Hall drift 
can transport magnetic field from the inner crust where the Ohmic decay is 
slow, to the outer crust where it proceeds rapidly. Unlike the Ohmic decay, 
the Hall drift is essentially a non-linear process. 

As we can see, the characteristic time scales of field decay strongly de- 
pend on the length scale. This brings us to the problem of magnetic field 
localization. If the magnetic field is distributed in the entire NS, includ- 
ing the core, then these time scales are long so that the magnetic field de- 
cay is insignificant on time scales about few million years. In the modern 
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models (lAguilera et all 120081 ; iPons et al.l . [20091) it is assumed that the main 



part of the magnetic field is localized in the crust of a NS. In this case 
r 0hm ~ 10 5 - 10 6 yrs and r Hall ~ (10 2 - lO^^o/lO^G)- 1 



yrs^ 



The evolution of magn etic field according to the model by lAguilera et al 



(120081 ); IPons et al.l (120091 ) is well approximated by the following equation: 

exp(-t/r hm) 



B(t) = B min + (B — B n 



(3) 



1 + T hm/THall{l ~ exp(-t / ' T hrn)) ' 

Here Bq is the initial magnetic field strength and B m i n = min[Bo/2 ,2 x 
10 13 ] G. It can be seen that an initially strong magnetic field evolves to 
the asymptotic value B = 2 x 10 13 G, which can be slightly different for 
different objects. It is probably so in the case of magnetar SGR 0418+5729, 
which was born wi th high magnetic field s trength and evolved to the value 
B < 7.5 x 10 12 G (IRea and Espositol . l201lh . For low-field objects the decay 
is not so pronounced and the field decays just by a factor of 2. Therefore 
according to the theory, there should be only a small number of objects with 
strong magnetic fields among populations of sources with ages ~ 10 6 yrs, for 
example among HMXBs. 

It is convenient to illustrate the evolution of the magnetic field on the 
P — P diagram. The results of evolution according to the eq. (|3]) for different 
T ohm and ThoXI are shown in four panels of FigfTJ In each panel evolutionary 
tracks are shown for different initial magnetic fields from 10 12 G to 10 16 G. 
The observed populations indicated in the plots are the following: pulsars 
(black dots), magnetars (open diamonds), five of the M7 (triangles), and 
RRATs (filled diamonds). According to the simple model of field decay 
(eq. [3]) sooner or later magnetic fields of initially strongly magnetized NSs 
converge to the asymptotic value B ~ 2 x 10 13 G, and clustering of sources 
can be visible. For PSRs, M7, or any other population this effect is not 
visible in the data. Therefore we can assume that this convergence should 
happen at larger ages, where no isolated NSs are observed now. 

The plot for a more realistic set of parameters is shown in the top left 
panel: Tohm = 10 6 yrs and THaii = 10 4 (£>o/10 15 G) _1 yrs. Therefore, we use 
these parameters for theoretical estimates in sections 5 and 6. At the other 
panels the convergence is too rapid to be consistent with the observational 
picture. 



1 Below we discuss the field values at the magnetic pole: B = B po i e — 2B equat0 r — 
2/i/R 3 , where R is the radius of a NS. For numerical estimates we use the value R = 10 6 cm. 
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Figure 1: P - P diagram of radio pulsars, magnetars and radio-quiet isolated NSs. The 
observed populations indicated in the plots are the following: pulsars (black dots), mag- 
netars (open diamonds), five of the M7 (triangles), and RRATs (filled diamonds). In each 
panel NSs evolutionary tracks are shown according to analytical approximation of the 
field decay model (Eq. 3) with different initial magnetic fields from 10 12 G to 10 16 G. 
In different panels different Tohm and THall are used. Top left: Tohm = 10 6 yrs, mall = 
lO^Bo/lO^G)" 1 yrs. Top right: r ohm = 10 7 yrs, T HaU = 10 2 {B Q /W*G)- X yrs. Bottom 
left: T 0hm = I0 6 yrs, T HaU = 10 3 (So/10 15 G)- x yrs. Bottom right: r ohm = I0 5 yrs, 
THaii = 10 3 (Bo/10 x5 G)- 1 yrs. Tracks are plotted for B = I0 12 , 10 13 , 10 14 , 10 15 , and 
10 16 G. In all cases NSs are treated as orthogonal rotators according to the magneto- 
dipole formula. Long-dashed lines correspond to ages I0 5 and 10 6 yrs. 



The convergence must occur beyond the death line. Then PSRs cannot 
be used to test this effect. The M7-like sources are too cold when the conver- 
gence appears. Magnetars are not active any more when the field reaches to 
the asymptotic value. We propose that determination of the magnetic fields 
of NSs in HMXBs can be used to see the effect of convergence on the time 
scale of few million years. We use Be/X-ray binaries to test predictions of 
theoretical models of field decay. 



3. Methods of magnetic field evaluation 



In this section we describe several methods to estimate the magnetic field 
of accreting NSs used in this paper. At first, we discuss the equilibrium pe- 
riod hypothesis in the case of wind and disc accretion. Then, we formulate 
the models based on the values of maximal observed rates of spin-down and 
spin- up of an accreting NS. After that we present the model based on the 
assumption of the equality of Alfven and corotation radii at the moment 
of appearence of pulsations. Finally, we summarize th r ee m ore sophist i- 
cat ed models: the m o del by llllarionov and Kompaneets (ll990T ). the model 
by iBisnovatvi-Koganl (ll99ll ). and a new approach based on the model by 



Shakura et all (1201 if ). 



3.1. The equilibrium period hypothesis 

In an accreting system a NS is a f fected by accelerating and decelerating 
torques (see, for instance, iLipunovl (119921 ) to which we refer for details of 
accretion physics basics). The spin- up of a NS is due to the orbital angular 
momentum brought by the accreted matter. In the case of disc accretion 
the accelerating torque can be represented as a flow of angular momentum 
transported from the last Keplerian orbit: 



K su , disc = M\J GMeRA- 
For wind accretion the following equation can be used: 



K 



su.wind 



Mi]VtR 2 G . 



(4) 



(5) 



Here M is the accretion rate, R A = (/z 2 /2M\/GM) 2/7 is the Alfven ra- 
dius, M is the NS mass, Q = 2n/P orb is its orbital frequency, and Rc = 
2GM / (vl rb + v 2 ind ) is the Bondi radius, v winc i is the stellar wind velocity, v or b 
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is the orbital velocity, e, 77 - some numerical coefficients, whose values depend 
on the properties of the acrretion flow. 

To determine the Bondi radius, Rg, we use the orbital velocity of a NS at 
periastron. For the ec centricity w e use the value e = 0.3 as the typical one 
for Be /X-ra y systems (|Reia. 120111). We tak e the stellar wind velocity v W i n d = 
v (r I ' R^fl 2 (IRaguzova and Lipunovl . Il998l ). where v is the wind velosity at 
the surface of a Be-star, r is the radius of a circumstellar outflowing disc, 
and R* is the Be-star radius. We use the following values: Vq = 10 km s _1 , 
r/R* = 10. 

The decelerating torque arises due to the magnetic field interaction with 
the accr etion flow. Thi s torque at the stage of accretion can be estimated as 
follows ( jLipunovl . Il992l ) : 



K 



sd 



~ kt w 



(6) 



Here \i is the mag netic moment of a NS and R co = (GMP 2 /A-n 2 ) 1 ^ is the 
corotation radius, P - the spin period of NS. k t is a numerical coefficient 
dependent on the model of interaction between the magnetoshere and sur- 
rounding matter. 

The rate of changes of the total angular momentum of a NS can be written 
down as follows: 



dluj 



My/GMeR A - k t ^ 2 /R, 
Mr]QR 2 G - k t fi 2 /R 



3 

co • 

.3 

co ■ 



for disc accretion 
for wind accretion 



(7) 



10 45 g cm 3 is the NS moment of inertia and oj 



2-k/P - is the 



where / 
spin frequency. 

The equilibrium hypothesis states that an accreting NS spends most of 
the time in the state with decelerating and accelerating moments nearly 
equalized. When accretion starts the system rapidly evolves towards the 
equilibrium state (see also section This state is characterized by a char- 
acteristic spin period, P eq . To calcu late it we use the following values of 



numerical coefficients ( Lipunov . 19921 ): e = 0.45, 77 = 1/4, k t = 1/3. For the 
case of disc accretion the equilibrium period then is the following: 



eq.disc 



2 15 / 14 vrA; t 1/2 e - 1 /V 6/7 M _3/7 (G'M)- 5 / 7 = 5A7fi 30 e / 7 M w 3/7 s. 
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For wind accretion: 



P eq . mnd = J ^P^v 2 (GM)~ z/2 M~ 1,2 ii = 47.29(P orb /10days) 1/ W 8 2 Mi6 1/ ^ 

(9) 

Here M\q is the accretion rate normalized to a typical value 10 16 gs 1 , /i3o is 
the magnetic moment normalized to a typical value 10 30 G cm 3 , and v$ = 
v/10 8 cms _1 is the velocity of the accreting matter. 

Assuming that the spin period of a NS is equal to its equilibrium period, 
the magnetic field B for disc accretion can be estimated as follows: 

B = 2- l ' A ix- 7 ' G kf ,12 t 7 l 2A P 7 "M l l 2 {GMfl & R-\ (10) 

For wind accretion: 



B = 2Jp-p- r l /2 v- 2 (GMf 2 M l ' 2 PR-\ (11) 
V k t n 

3.2. Estimate based on the maximum spin-down rate 



This method is described, for instance, in iLipunovl (119921 ) . 

Observations show that NSs in binaries exhibit episodes of deceleration 
and acceleration. At the moment of the maximum spin-down the decelerating 
torque is much larger than the accelerating one. In this case we can neglect 
the first term in eq. ((Tj) and obtain the following: 

dlui , a 2 

-h-^r- (12) 



dt -R 3 



CO 



Using the observed values of the maximum spin-down rate the magnetic 
field of a NS can be estimated as follows: 

2 IPGM] 

This estimate should be normally considered as a lower limit, since we cannot 
be sure that no accelerating torque exists at that moment. 

The data to distinguish episodes of the maximum spin- down and to deter- 



mine the value of P is taken from graphs in lGalache et al.l ( 120081 ). As sources 
are transient and accretion occurs mainly close to the periastron passage, we 
use only data points from the same episode of accretion to determine the 
maximum spin-down rate. 
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3.3. Estimate based on the maximum spin-up rate 

The description of this approach also can be found in iLipunovl (ll992f ). 

At the moment of the maximum spin-up the accelerating torque is much 
larger than the decelerating one. In this case we can neglect the second term 
in eq. (j7j) and obtain the following: 

^ = M{GMeR A ) l/2 ■ (14) 

This equation is written for the case of disc accretion. Using the observed val- 
ues of the maximum spin-up rate, the magnetic field of NS can be estimated 
as follows: 

B ^ ( /P ) 7/2 (15 ) 

e 7/4 R3 P 7M3(GM) S / 2 ' 

This estimate should be normally considered as a lower limit to the value 
of the magnetic field, since we cannot be sure that no decelerating torque 
exists at that moment. 



We use graphs in iGalache et al.l (120081 ) to find episodes of the maximum 



spin-up and to determine the value of P as described in the previous subsec- 
tion. 

3. 4- Onset of accretion and the condition Ra = R co 

In our work we consider only Be/X-ray binaries. Many systems of this 
type have large eccentricities. During its orbital motion a NS can switch from 
the propeller to the accretor phase. At this moment the condition Ra = R co 
is roughly satisfied: 

jLjf-f™*^ 1 *. (16) 



^My/GMJ V 47r2 

We assume that at this moment pulsations start being registered. Using 
spin period and luminosity we can estimate the magnetic field using the eq. 
ffTB]) as follows: 

B = 2-^\- 7 ^R~ 3 (M) 1 / 2 (GMf 6 P 7 / 6 . (17) 
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3.5. The model by Illarionov and Kompaneets 

Above we described the well -known approa ches developed by different 
authors in 70s-80s (see a review in Lipunov ( 19921 )). No w we come to more so- 



phisti cated models, starting with the model proposed by llllarionov and Kompaneets 



(1990) 



This is a model with a more efficient spin-down mechanism of an accreting 
NS. The spin-down is a result of an efficient angular momentum transfer from 
the rotating magnetosphere of an accreting object to the outflowing stream 
of magnetized matter. The outflow is formed within a limited solid angle, 
and the outflow mass rate is less than the accretion rate. Frozen magnetic 
fields in the outflow may provide additional angular momentum losses. 

We assume, that at the moment of the maximum spin-down the decel- 
erating torque is much larger than the accelerating one. This decelerating 
torque is estimated according to the model by Illarionov and Kompaneets. 

The angular m omentum losses in a solid angle y of an outflow can be 



written as follows (jlllarionov and Kompaneets! . Il990h : 



J^L = -k^MuR A . (18) 

at Z7T 

We use this approach to obtain the decelerating torque in a small solid angle 
X ~ 1- For the rest of the NS magnetosphere we assume that the spin-down 
mechanism described above (eq. fT2l) is valid. Then the equation for the 
spin evolution of a NS can be written in the following way, neglecting the 
accelerating moment: 

^ = -k—MuR 2 A - k t ^' X y . (19) 
dt 2vr A 4nR 3 co v ' 

To our knowledge, there is no analytical solution for this equation for mag- 
netic field strength. Therefore we solve it numerically using bissection method. 
For numerical estimates below we use the following parameter values: x = 1, 
k = 2/3, and k t = 1/3. We estimate a NS magnetic field fr om this equa- 



tion u sing measured P and P which are taken from graphs in iGalache et al. 



((20081). 



3.6. The model by Bisnovatyi-Kogan 



In the model by iBisnovatyi-Koganl (119911 ) the author proposed to evaluate 



spin periods of NSs as: P = \JP\P%. Here P\ is the equilibrium period of a 
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NS. The decelerating torque calculated according to the model by Illarionov 
and Kompaneets (eq. [TSJ) is compensated by an accelerating torque (see eq. 



orb 



P( V 2GM 



Here we use the Alfven radius in the form R A = (B 2 R 6 /(8MVGM)) 2/7 , 
where B is the field at the magnetic pole. The numerical coefficient k is 
assumed to be 2/3. 

The equilibrium period in this case can be calculated as follows: 

P = ^g^gg m) 
1 2 11 /4 V ^(M)V2(G'M) 3 / 2 

Another characteristic period, P2, is estimated in assumption that the 
angular velocity of the NS is equal to the angular velocity of the infalling 
matter on the distance of half Alfven radius: 

This equation differs from the equation given in (iBisnovatyi-Kogan . Il99l|) by 
a factor of two. 

From equation (1221) the period P2 can be estimated as: 

v^P orb B^R^ 



8 4 / 7 16M 4 /7(GM) 16 / 7 
Therefore, the spin period of a NS can be calculated as follows: 

P= /p^ = 0,37 orb ^— . (24) 

Using relation f )24|) we can estimate the magneic field as follows: 



.P u / 15 My 2 {GM) 



53/30 



B = 2.95 . (25) 

v^P 7 J b 10 R* 
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3. 7. Hybrid model 

Recently, a ne w detailed mod e l for the case of quasi-spherical accretion 
was developed by IShakura et al.l ( 1201 lh . In this model a settling subsonic 
accretion proceeds through a hot shell formed around a NS magnetosphere. 
The turbulent stresses are capable of carrying the angular momentum out- 
wards through the shell. 

In this way it is possible to estimate the magnetic field of a NS as follows: 



B 12 = 0.24 



-12/33 \ n /12 



l-z/Z 



P/lOOs V 1/12 ^i/a 11/3 
P or6 /10daysJ Ml6 V8 ■ 



(26) 



12 



Here we assume that the spin period is equal to the equilibrium value and 
(B/10 12 ) G. The expr ession in the b r acket s (5 2 ^" 12 / 33 )/(l - z/Z) 



B 

is approximately equal to 1 (IShakura et al. 
numerical coefficients. 



201lh . 5,m,£,z, and Z - are 



The model developed by IShakura et al.l (120111 ) is applicable to long pe- 
riod NSs with quasi-spherical accretion. Accretion in Be/X-ray binaries can 
proceed with or without formation of an accretion disc around a compact 
object. It is important to distinguish objects where an accretion disc can be 
formed. Th ere is a well-known condition for accretion disc formation (see, 
for instance, Illarionov and Sunyaev ( 1975h ) : 



r]QR 2 G > (GMR A 



,V2 



(27) 



The values of numerical coefficient r\ depend on the structure of stellar wind, 
presence of inhomogeneities and other properties of accreting matter. For 
large values of rj ~ 1 an accretion disc is formed for all of the systems 
under discussion. iDavies and Pringld ( 119801 ) showed that the coefficient can 
be much smaller than 1. In this case a n accretion disc may fail t o form for 
systems with large values of spin period (IDavies and Pringld . 119801 ). We take 
into account this condition for the model by Shakura et al. We apply this 
model only to objects with P > 50 s. For short-period objects (P < 50 s) 
we make calculations with the hypothesis of the equilibrium period using 
standard formulae for disc accretion. So, we call this model "the hybrid". 



4. Observational data 



For our estimates we use data (IGalache et al.l . 120081 ) on 40 Be/X-ray bi- 
naries in the SMC observed by the RXTE satellite. These data include: the 
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spin period and count rate (in counts PCU x s 1 ) for the pulsating compo- 
nent of emission. We convert counts to flux using WebPIMMdl Besides the 
standard RXTE PCA parameters we need to specify the Galactic hydrogen 
column density (uh) towards the SMC and spectral parameters. 

For the Galactic hydrogen colum n density in the direc tion of the SMC 



we take the value n H = 4 • 10 cm (IKuulkers et al.l . 120071 ). Our results are 



not very sensitive to this parameter. Even a significant variation of n# by 
a factor of a few results in the magnetic field variation on the level of a few 
percent. 

We assume that spectra can be well approximated by a single power-law. 
So, the only spectral parameter is the photon i ndex. The photon index of 
all sources is taken to be equal to 1 (jLipunovl . Il992f ) . The magnetic field 
estimates show only modest sensitivity to this parameter. When we change 
a photon index by a factor of 1.5, the magnetic field changes by ~ 10%. 

For our magnetic field estimates we need to know the value of the ac- 
cretion rate, M. We derive it from the X-ray flux value. Formally, we have 
to use the bolometric flux, however we use a simplified approach. From the 
measured count rate we derive the flux in some output energy range using 
WebPIMMS, assuming a single power- law spectrum. To estimate the out- 
put energy range we co nsider the data for different Be /X -ray systems in the 
Galaxy from the paper (jLutovinov and Tsygankovl . 120081 ) . In this work it is 
shown that most of the objects have power-law spectra up to 25 keV, and 
we accept this value f or our study. The low er limit - 3 keV - for the output 
energy is taken from ( Galache et all 120081 ) . So, we have the output energy 
range 3-25 keV to est i mate the luminosity used in M calculations. 

In (1 Galache et all 120081 ) only the data on pulsed flux are given. For our 
estimates we need to make assumptions about the total flux, as finally we 
need to derive the value of the mass accretion rate. Obviously, the ratio be- 
tween pulsed and total fluxes may be different for different sources. However, 
after a statistical study we decide to use a single value for the whole set of 
sources, as no details are available for most of the objects in the sample. We 
assume that the total flux is 5 times larger than the pulsed flux (Lutovinov, 
private communication) . 



2 http://heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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5. Results 



Our main results consist of estimates of NSs magnetic field in Be/X-ray 
binaries obtained by several different methods described in Sec. 3. They are 
presented in Figs. |2] and [3J Estimates for different models are shown with 
different symbols. In Fig. [5] estimates of magnetic field with the assump- 
tion of the equilibrium period hypothesis are plotted as filled circles, and as 
open circles in cases of wind and disc accretion, respectively. Estimates for 
the hypothesis of the maximum spin-down and spin-up are drawn as down 
triangles and up triangles, correspondingly. In Fig. [3] with filled triangles 
we show estimates of magnetic fields with the assumption of the onset of 
pulsations when Ra = R co , and for the model by Illarionov and Kompaneets 
- as diamonds. Finally, for the model by Bisnovatyi-Kogan and the hybrid 
model our results are given as crosses and asterisks correspondingly. 

It is expected that the spin periods are on average larger for NSs with 
larger magnetic fields. This is due to the fact that for a given spin period 
for larger magnetic field the decelerating torque is much larger, see eq. fl7J), 
but the accelerating one remains unchanged (for quasi-spherical accretion), 
or just slightly increases (if an accretion disc is formed). So, the equilibrium 
is reached at longer spin periods for larger fields. The accelerating moment 
linearly depends on the accretion rate; in the case of wind accretion without 
a disc it also depends on the properties of the binary. This normally should 
result in a scatter in the relation between spin periods and magnetic fields. 

In Figs. |2] and [3] we see that for some models the correlation between 
spin period and magnetic field of NS is very tight. This can be a result of 
a selection effect. In our case pulsations are registered only above some flux 
value. All sources are nearly at the same distance. Then we immediately 
obtain that in our approach there is a critical value of M, below which 
pulsations are undetectable. This effect is especially important for the model 
of equilibrium period, for the hybrid model, and for the onset of accretion 
when Ra = R co - Therefore, this effect makes the correlation B — P very 
tight reducing the scatter. 

There is a well-known problem that for long-period objects traditional 
methods (equi librium period, Ra = R co , and the maximum spin-down, see 



Lipunov! (119921 )) often produce too large field value estimates. This is clearly 



visible in Figj2j On the other hand, the model of the maximum spin-up 
shows more realistic results. Still, this m ethod strongly depen ds on mea- 
sured parameters. Using the data given in iGalache et al.l (120081 ) we cannot 
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Figure 2: Magnetic field values estimated using different methods: filled circles - the 
equilibrium period hypothesis in the case of wind accretion, open circles - the equilibrium 
period hypothesis in the case of disc accretion, down triangles - the maximum spin-down, 
up triangles - the maximum spin- up. 
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Figure 3: Magnetic field values estimated using different methods (continued): filled tri- 
angles - onset of pulsations when Ra = R C o, asterisks - the hybrid model, diamonds - the 
model by Illarionov and Kompaneets, and crosses - the model by Bisnovatyi-Kogan. 
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measure period variations and corresponding time intervals with high preci- 
sion. The spin periods are presented in the plots with large error bars, and 
time intervals corresponding to episodes of the maximum spin-up are very 
short. Therefore the values of P are obtained with significant uncertainties, 
so estimates obtained with this method are not very certain. However, we 
present results based on the maximum spin-up for long-period sources for 
which estimates seem to be reliable. 

Estimates based on the model by Bisnovatyi-Kogan are more reliable 
than those based on the maximum spin-up, and do not contain very high 
values > 10 14 G. Therefore, there are many NSs with B > 3 x 10 13 G in 
contradiction with theoretical expectations. 

Finally, the hybrid model is able to reproduce most of the predicted deficit 
of long-period high magnetic field objects. Thus, we consider this model to 
be the best among studied and discuss results based on it in the next section 
in more details. 



6. Discussion 

In Fig. 4 we show magnetic field distributions. Here the theoretical 
distribution for the initial polar magnetic field -Bo is generated randomly for 
10 5 NSs from a log-normal distributon with the mean value < lo g Bq/[G] >= 



13.25 and standard deviation 0i og B o = 0.6 (IPopov et all 120101 ) . Note, that 
formally this initial field distribution is not absolutely self-consistent with 
our model of field decay, because we use parameters slightl y different from 



those used in numerical calculations by ( IPopov et all l2010f ). However, this 
imperfection cannot significantly influence our conclusions, because the main 
assumption is related to the asymptotic field value. 

The evolved fields are calculated according to eq. (j3J) using the values: 
THaii = 10 4 (-B /10 15 G) -1 yrs and T 0hm = 10 6 yrs. The NS ages are chosen 
uniformly in the interval [0,t maa! ], where t max = 10 7 yrs is the maximum life 
time of Be/X-ray binaries. Magnetic field estimates obtained in this paper 
are shown in four different panels in Fig. HI In the panel A with a solid line 
we show the magnetic field distribution based on the assumption of the onset 
of pulsations when Ra = R co - The magnetic field distribution according to 
the model by Bisnovatyi-Kogan is plotted in the panel B. Distributions for 
the models by Illarionov, Kompaneets and the hybrid model are shown in the 
panels C and D, correspondingly. Theoretical models for the initial and the 
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Figure 4: Magnetic field distributions for different models (solid lines). Panel A - onset 
of pulsations when Ra — R C o- Panel B - the model by Bisnovatyi-Kogan. Panel C - the 
model by Illarionov and Kompaneets. Panel D - the hybrid model. Theoretical models 
for the initial and evolved magnetic fields are shown in all the four panels by dashed and 
dash-dotted lines, respectively. All distributions are normalized to the total number of 
observed objects. Error bars are given according to the Poissonian values. 

evolved magnetic fields are shown in all four panels by dotted and dot-dashed 
lines, correspondingly. 

To compare the theoretical predictions with results based on the hybrid 
model we use the Kolmogorov-Smirnov test. According to this test, the 
hypothesis of the same distribution of the theoretical and estimated values 
cannot be rejected at 80% significance level. However, there is a small number 
of objects for the hybrid model distribution in the bin 10 13 < B < 10 13 ' 5 G 
(Fig. HP). Discrepancies can be related to uncertainty in evolution of the 
magnetic field for t > 10 6 yrs. Probably, the magnetic field decay continues 
but more slowly than before. Likewise, different highly magnetized objects 
can evolve to different values of the asymptotic magnetic field (see Sec. 2). 
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In this study we assumed that NSs in binaries have the same initial mag- 
netic distribution as isolated compact objects. This is not so obvious in 
the case of magnetars, as in the starda rd scenario their magne t ic fiel ds are 
generated via the dynamo mechanism (IDuncan and Thompson! . Il992l ). The 
dynamo mechanism is effective in the case of very rapid rotation of a pro- 
toNS. Isolated massive stars hardly can have rapidly rotating cores. Rapid 
rotation can be reached in close binary systems due to accretion or/and tidal 
synchronization. On the other hand, all known magnetars are isolated ob- 
jects, and authors studied evolutionary channels in binary systems which 
result is formation of isolated compact objects with rapid initial rotation 



(IBogomazov and Popovl . 120091 ; iPopov and Prokhorovl . 120061 ). In these mod- 



els objects with extremely large field tend to form after a binary is destroyed 
to explain that all known magnetars are isolated. 

To test the possibility of absence of extreme magnetars in binary systems 
for the case of the hybrid model (which is the most successful in explaining the 
properties of the population under study), we plot an additional histogram 
(Fig. \5§ for a modified initial magnetic field distribution: no NSs with B > 
10 14 G are formed. As we see, estimates made on the basis of observational 
data are even in better agreement with the prediction for the evolved field 
in the model by iPons et al.l (120091 ) with the cutoff at Bq = 10 14 G, as the 
discrepancy in the bin 10 13 < B < 10 13,5 G is decreased. However, we note 
that our calculations are b ased on an approximation to the actual detailed 
model by lPons et al.l (120091 ). and for the specific choice of parameters, which 
can vary (rHaihTohm, B min ). We do not overestimate the accuracy of our 
results, and so do not make further speculations. 

An additional consisten cy check for m agnetic field estimates can be made 
using the Corbet diagram (jCorbetl . Il984l ). We apply this test to the hybrid 
model. This diagram relates spin periods and orbital periods of NSs. We 
simulate this distribution. Orbital periods of NSs are chosen uniformly in 
the interval [10, 1000] days. Using an orbital period and evolved magnetic 
field we can estimate the spin period of a NS. We calculate the evolved 
magnetic field as discussed above. These values of spin and orbital peri ods 
are plotted in Fig. 6 as dots. We plot the observational data (IReigl . 1201 if ) as 
diamonds in this figure. One can see good agreement between theoretical and 
observational data. According to the non-parametric Kolmogorov-Smirnov 
test, the hypothesis that both subsets arise from the same distribution cannot 
be rejected at 77% significance level. 

The methods that we use to estimate magnetic fields have limitations, 
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Figure 5: Magnetic field distributions. The distribution based on observational data after 
application of the hybrid model is shown with a solid line. Poissonian error bars are 
indicated. Results for the theoretical model for evolved magnetics fields is shown with a 
dash-dotted line. Finally, a modification of the theoretical model for which a cut-off at 
10 14 G for initial fields was used is shown with a dashed histogram. All distributions are 
normalized to the total number of observed objects. 
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especially the most simple ones. For instance, the hypothesis of the equi- 
librium period may not be valid for some systems. In our estimates we use 
M obtained for the moment when pulsations are detected for the first time. 
But for NSs in systems with high eccentricity which spend much time at 
the propeller stage, it can take some time to reach the equilibrium. So, 
their observed periods are larger than the equilibrium value. This results in 
underestimating of the magnetic field. 

Likewise, we should consider that both disc and wind accretion estimates 
(see, Fig. [2]) cannot be correct at the same time because just one regime 
is realized in any given source. On average, estimates for disc accretion are 
more realistic for shorter spin periods. Still, we think that estimates based 
on old models are not realistic in general, as they produce too large field 
values for long peri od sources. 

In the model by IShakura et al.l (1201 if ) there is a critical value of M. For 
accretion rates above this M cr an envelope around a compact object can cool 
rapidly. This results in an enhanced accretion rate, and no equilibrium is 
possible in this regime. This critical value depends on several parameters 
and can be different in different systems, but roughly M cr < 10 17 g/s for 
magnetic fields > 10 13 G. Marginally, this condition is fulfiled for our systems. 
We can slightly over estimate the accretion rate, however, if the model by 



Shakura et al.l (120111 ) is applicable, then the magnetic field estimate depends 



on the accretion rate only mildly. 

It is instructive to compare the results of indirect magnetic field esti- 
mates with direct magnetic field measurements through cyclotron lines, as 
it was made for the long-per iod object GX 301-2. Fo r this object, magnetic 
field strengh is estimated by lLa Barbera et al.l ( 120051 ) as ~ 5 x 10 12 G. Th is 



value is consistent with the magnetic field estimate made by IShakura et al 



(120111) using their mode l. On the other hand, the classical models used by 



Doroshenko et al.l ( 120101 ) provide the magnetic field estimate nearly two or- 



ders of magnitude larger than it is given by the cyclotron line measurements. 
However, the cyclotron line can be formed high above the NS surface that 
allows the measured magnetic field str e ngth to be considerably weaker than 



that on the surface (IDoroshenko et al.l . 120101 ) . 



Finally, we comment on a possible evolution of the angle between spin 
and magnetic axis. For radio pulsars it is possible to d escribe the a ngular 
momentum evolution using the magneto-dipole formula (jPacinil . 119671 ): 

dlu lB 2 u 3 R 6 „ 

15 — sin 2 /3, (28) 



dt 
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Figure 6: Corbet diagram for the observed data (open diamonds) and for the hybrid model 
(dots). 
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Here (3 is the angle between spin and magnetic axis. From this formula one 
can estimate the effective magnetic field -B e // — -Bsin/3. So, the effect of 
the variation of the effective magnetic field can be related not only to the 
magnetic field decay, but also to evolution of the angle between magnetic 
and spin axis. Therefore, it is difficult to distinguish between these effects. 
However, there are different methods to estimate magnetic field of magne- 
tars and to demonstrate that this is really the field decay, not the variation 
of the angle 0. Note, that our methods are only weakly sensitive to the 
angle between spin and magnetic axis, and we estimate the absolute value 
of magnetic field. Still, using Be/X-ray systems we can try to study evolu- 
tion of the angle between these axis. According to two most popular models 
(magneto-dipole and current losses), the angle evolves very fast on the time 
scale of spin period evolution. The angle (3 evolves to the value = 0° in 
the ma gneto-dipole model a nd to the value (3 = 90° in the current losses 
model (jGurevich et al.l . Il993l ). But this effect is not observed for PSRs. It is 
possible to use Be/X-ray binaries to study evolution of the angle on the time 
scale up to 10 7 yrs jointly with the magnetic field decay. We plan to use this 
approach in future studies. 



7. Summary 

In our work we used different models to estimates magnetic fields of NSs 
in Be/X-ray binaries in the SMC. Most of them are shown to o yerestimate the 



magne tic field, but the hybrid model based on the approach by lShakura et al. 



(120 111 ) gives good agreement with the prediction of the theoretical model of 
magnetic field decay on a time scale up to ~ 10 7 yrs. 



Acknowledgements 

We thank N.I. Shakura and K.A. Postnov for the opportunity to use their 
model before publication, and for consultations. We acknowledge A. A. Lu- 
tovinov and J. A. Pons for interesting discussions. We thank P.K. Abolmasov 
and E.A. Vasiliev for useful remarks. Also we are grateful to the anonymous 
referees for useful comments. Sergei Popov thanks University of Padova for 
support and hospitality during his visit. This work was supported by the 
Russian Foundation for Basic Research grants 10-02-00599, 12-02-00186 and 
by the Federal program for research staff (02.740.11.0575). 



23 



References 

Aguilera, D.N., Pons, J. A., Miralles, J.A., 2008. The Impact of Magnetic 
Field on the Thermal Evolution of Neutron Stars. ApJL 673, L167-L170. 
larXiv: 0712. 13531 

Bhattacharyya, S., 2009. X-ray views of neutron star low-mass X-ray binaries. 
Current Science 97, 804-820. larXiv: 1002.44801 

Bisnovatyi-Kogan, G.S., 1991. Rotational equilibrium of long-periodic X-ray 
pulsars. A&A 245, 528-530. 

Bogomazov, A.I., Popov, S.B., 2009. Magnetars, gamma-ray bursts, and very 
close binaries. Astronomy Reports 53, 325-333. larXiv : 0905 . 3238 . 

Corbet, R.H.D., 1984. Be/neutron star binaries - A relationship between 
orbital period and neutron star spin period. A&A 141, 91-93. 

Davies, R.E., Pringle, J.E., 1980. On accretion from an inhomogeneous 
medium. MNRAS 191, 599-604. 

de Luca, A., 2008. Central Compact Objects in Supernova Remnants, in: 
C. Bassa, Z. Wang, A. Cumming, & V. M. Kaspi (Ed.), 40 Years of Pulsars: 
Millisecond Pulsars, Magnetars and More, pp. 31 1-319. larXiv : 0712 . 22091 

Doroshenko, V., Santangelo, A., Suleimanov, V., Kreykenbohm, I., Staubert, 
R., Ferrigno, C, Klochkov, D., 2010. Is there a highly magnetized neutron 
star in GX 301-2? A&A 515, A10. 10907.38441 

Duncan, R.C., Thompson, C, 1992. Formation of very strongly magnetized 
neutron stars - Implications for gamma-ray bursts. ApJL 392, L9-L13. 

Galache, J.L., Corbet, R.H.D., Coe, M.J., Laycock, S., Schurch, M.P.E., 

Markwardt, C, Marshall, F.E., Lochner, J., 2008. A Long Look at the 

Be/X-Ray Binaries of the Small Magellanic Cloud. ApJ Supp. 177, 189- 
215. arXiv: 0802 .21181 

Goldreich, P., Reisenegger, A., 1992. Magnetic field decay in isolated neutron 
stars. ApJ 395, 250-258. 



24 



Gurevich, A., Beskin, V., Istomin, Y., 1993. Physics of the Pulsar Magneto- 
sphere, pp. 432. ISBN 0521417465. Cambridge, UK: Cambridge University 
Press, August 1993. 

Haberl, F., 2007. The magnificent seven: magnetic fields and surface 
temperature distributions. Astrophys. & Space Sciences 308, 181-190. 
larXiv: 06090661 

Illarionov, A.F., Kompaneets, D.A., 1990. A Spin-Down Mechanism for Ac- 
creting Neutron Stars. MNRAS 247, 219. 

Illarionov, A.F., Sunyaev, R.A., 1975. Why the Number of Galactic X-ray 
Stars Is so Small? A&A 39, 185. 

Kaplan, D.L., 2008. Nearby, Thermally Emitting Neutron Stars, in: 
C. Bassa, Z. Wang, A. Cumming, & V. M. Kaspi (Ed.), 40 Years of Pulsars: 
Millisecond Pulsars, Magnetars and More, pp. 331-339. 

Keane, E.F., 2010. Transient Radio Neutron Stars larXiv : 1008 . 36931 

Kuulkers, E., Shaw, S.E., Paizis, A., Chenevez, J., Brandt, S., Courvoisier, 
T., Domingo, A., Ebisawa, K., Kretschmar, P., Markwardt, C.B., Mowlavi, 
N., Oosterbroek, T., Orr, A., Rfsquez, D., Sanchez- Fernandez, C, Wij- 
nands, R., 2007. The INTEGRAL Galactic bulge monitoring program: 
the first 1.5 years. A&A 466, 595-618. larXiv: 070 12441 

La Barbera, A., Segreto, A., Santangelo, A., Kreykenbohm, I., Orlandini, 
M., 2005. A study of an orbital cycle of GX 301-2 observed by BeppoSAX. 
A&A 438, 617-632. 

Lipunov, V.M., 1992. Astrophysics of Neutron Stars, XIII, 322 pp. 108 figs. 
Springer- Verlag Berlin Heidelberg New York. 

Liu, Q.Z., van Paradijs, J., van den Heuvel, E.P.J. , 2006. Catalogue of high- 
mass X-ray binaries in the Galaxy (4th edition). A&A 455, 1165-1168. 
larXiv: 0707. 05491 

Lutovinov, A., Tsygankov, S., 2008. X-ray pulsars through the eyes of INTE- 
GRAL, in: M. Axelsson (Ed.), American Institute of Physics Conference 
Series, pp. 191-202. arXiv: 0808. 20341 



25 



Manchester, R.N., Hobbs, G.B., Teoh, A., Hobbs, M., 2005. The Aus- 
tralia Telescope National Facility Pulsar Catalogue. AJ 129, 1993-2006. 
larXiv: 04126411 

Mereghetti, S., 2008. The strongest cosmic magnets: soft gamma- 
ray repeaters and anomalous X-ray pulsars. A&ARv 15, 225-287. 
larXiv: 0804. 02501 

Negueruela, I., 2010. Stellar Wind Accretion in High-Mass X-Ray Binaries, 
in: J. Marti, P. L. Luque-Escamilla, & J. A. Combi (Ed.), High Energy 
Phenomena in Massive Stars, p. 57. arXiv T0907 . 2883 , 

Pacini, F., 1967. Energy Emission from a Neutron Star. Nature 216, 567-568. 

Pavlov, G.G., Sanwal, D., Teter, M.A., 2004. Central Compact Objects 
in Supernova Remnants, in: F. Camilo & B. M. Gaensler (Ed.), Young 
Neutron Stars and Their Environments, p. 239. arXi v: 0311 526, 

Pons, J. A., Geppert, U., 2007. Magnetic field dissipation in neutron star 
crusts: from magnetars to isolated neutron stars. A&A 470, 303-315. 
larXiv: 07032671 

Pons, J. A., Miralles, J. A., Geppert, U., 2009. Magneto-thermal evolution of 
neutron stars. A&A 496, 207-216. larXiv: 0812 .30181 

Popov, S.B., Pons, J.A., Miralles, J.A., Boldin, PA., Posselt, B., 2010. Popu- 
lation synthesis studies of isolated neutron stars with magnetic field decay. 
MNRAS 401, 2675-2686. larXiv: 0910. 21901 

Popov, S.B., Prokhorov, M.E., 2006. Progenitors with enhanced rotation 
and the origin of magnetars. MNRAS 367, 732-736. larXiv : 05054061 

Postnov, K.A., Yungelson, L.R., 2006. The Evolution of Compact Binary 
Star Systems. Living Reviews in Relativity 9, 6. arXiv: 0701059, 

Raguzova, N.V., Lipunov, V.M., 1998. High-eccentric X-ray binaries: evolu- 
tion, wind rose effect, accretor-propeller luminosity gap. A&A 340, 85-102. 

Raguzova, N.V., Popov, S.B., 2005. Be X-ray binaries and candidates. As- 
tronomical and Astrophysical Transactions 24, 151-185. arXiv: 0505275. 



26 



Rea, N., Esposito, P., 2011. Magnetar outbursts: an observational review, 
in: D. F. Torres & N. Rea (Ed.), High- Energy Emission from Pulsars and 
their Systems, p. 247. 11101.44721 

Reig, P., 2011. Be/X-ray binaries. Ap&SS 332, 1-29. larXiv: 1101 .50361 

Shakura, N., Postnov, K., Kochetkova, A., Hjalmarsdotter, L., 2011. Theory 
of quasi-spherical accretion in X-ray pulsars, accepted for publication in 
MNRAS I arXiv: 1110.37011 



27 



